HfC^: a new direction for intrinsic defect driven ferromagnetism 
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In view of the recent experimental reports of unexpected ferromagnetism in Hf02 thin films Q], we 
carried out first principles investigations looking for magnetic order possibly brought about by the 
presence of small concentrations of intrinsic point defects. Ab initio electronic structure calculations 
using density functional theory (DFT) show that isolated cation vacancy sites in Hf02 lead to the 
formation of high spin defect states. Furthermore these appear to be ferromagnetically coupled with 
a rather short range magnetic interaction, resulting in a ferromagnetic ground state for the whole 
system. More interestingly, the occurrence of these high spin states and ferromagnetism is in the 
low symmetry monoclinic phase of Hf02- This is radically different from other systems previously 
known to exhibit point defect ferromagnetism, warranting a closer look at the phenomenon. 

PACS numbers: 75.10.-b,75.50.-y,75.70.-i 



Pure Hafnium oxide (Hafnia, Hf02) is a wide band gap 
insulator with a high dielectric constant and presents no 
evidence of any magnetic order in the ground state. It 
has been heavily studied in recent times within a first- 
principles context 0, Q along with several other high-k 
dielectrics mainly because of its potential for substitut- 
ing SiC>2 as a gate dielectric in microelectronic devices. 
First-principles investigations of the physics of defects in 
Hafnia [4j,|5j also reflect this trend. 

Given the intrinsic non-magnetic nature of HfC>2, the 
discovery of ferromagnetic order in thin films which are 
pure except for the possible presence of intrinsic defects 
[l| comes as a surprise. Since the nominal valence of Hf in 
Hf02 is 4+, which leaves Hf atoms with an empty d-shell, 
the phenomenon was initially termed dP magnetism. It 
was suggested that the magnetism probably arose from 
partially filled d-orbitals derived from Hafnium atoms co- 
ordinating Oxygen Vacancy sites (Vo). In this letter we 
present the results of our first-principles band structure 
calculations which clearly show that the observed ferro- 
magnetism is most likely due to the presence of cation 
(Hf) vacancy sites Vgf. These form high spin states de- 
rived mainly from the Oxygen p-orbitals co-ordinating 
the Vkf- 

Intrinsic point defect driven ferromagnetism in other- 
wise non magnetic compounds has been previously stud- 
ied using first principles methods, in several systems. 
Notable cases include CaB 6 @], CaO Q and SiC @. 
The main characteristic of the vacancies in all these sys- 
tems is the high symmetry, either octahedral or tetra- 
hedral around the vacancy site. This invariably leads 
to a highly degenerate single particle spectrum, which 
may then present high spin states. In CaO for instance 
Coulomb repulsion stabilizes the two holes occupying the 
degenerate molecular orbital associated to Vc a in & triplet 
ground state . Similarly Zywietz et. al showed that for 
a Si vacancy in cubic SiC, only the occurrence of a mag- 
netic Jahn- Teller distortion stabilizes the spin singlet rel- 



ative to the triplet state otherwise expected from the Td 
symmetry and the degenerate single particle spectrum. 
Hafnium vacancies in monoclinic Hf02 however, are rad- 
ically set apart by the complete lack of symmetry around 
the vacancy site. In a defect molecular model, the sin- 
gle particle spectrum is completely non-degenerate and 
yet, as our DFT calculations show, they have magnetic 
ground states. This suggests that our understanding of 
defect induced ferromagnetism is at least incomplete and 
the phenomenon deserves a closer look. 

Our DFT calculations are carried out in the local spin 
density approximation (LSDA) using the Ceperly-Alder 
exchange correlation potential in the numerical im- 
plementation contained in the code siesta 10]. Test 
calculations using the generalized gradient approxima- 
tion (GGA), yield qualitatively similar results. We use 
the numerical localized atomic orbital basis set imple- 
mented in siesta including polarized orbitals with an 
energy shift of 0.01 eV (Toj. The basis set consists of O 
2s, 2p states and Hf 6s, 5d and 6p states. A set of two 
zetas and a polarized orbital are used for every shell ex- 
cept the 5d for which two zetas suffice, fc-point sampling 
is done by using a fixed fc-grid cutoff of 15.0 A which is 
equivalent to a 6x6x6 Mokhorst-Pack mesh for the twelve 
atom monoclinic unit cell. For 2x2x2 supercell contain- 
ing 96 atoms the same cut-off is equivalent to a 3x3x3 
Monkhorst-Pack mesh. The energy cutoff defining the 
equivalent planewave cutoff for the numerical grid is set 
at a value of 250.0 Ryd. In all cases, the systems stud- 
ied are relaxed until all the forces are smaller than 0.05 
eV/A. 

In order to make contact with the d° magnetism model 
we first look at the case of Oxygen vacancies. In the 
monoclinic structure there are two possibilities: i) Vq 
co-ordinated by three Hf atoms (V03), and ii) Vo co- 
ordinated by four Hf atoms (V04). These are shown in 
Fig. n A description of the various defect energetics is 
outside the scope of our work (see reference and here 
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we focus only on the electronic and magnetic structure 
of the defects. 




FIG. 1: The geometries around Oxygen vacancies V03, V04 
and a Hafnium Vacancy VHf are shown. The central atoms 
are included for clarity. The V04 site forms an imperfect 
tetrahedron while the V03 site forms an almost planar trigo- 
nal structure with the co-ordinating Hf atoms. 

An Oxygen vacancy in both the three and four fold 
co-ordinated case leads to the formation of a set of im- 
purity levels, with a low lying level in the HfC>2 band gap 
(V03 and V04) and higher lying levels just below the 
conduction band (V03* and V04*, see figure EJ). These 
impurity levels are formed from the dangling <i-orbitals 
of Hf atoms co-ordinating the vacancy site. The level 
in the band gap is filled by two electrons which in the 
perfect crystal, would have populated O p-orbitals. In 
this sense Vo is an n-type defect as it results in two 
electrons occupying conduction band derived cation or- 
bitals. Also from the density of states (DOS) it follows 
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FIG. 2: Spin resolved total DOS for Oxygen deficient Hf0 2 . 
The defect states are labeled according to figure Q V04* and 
V03* denote the empty states at the bottom of the conduc- 
tance band. The vertical line indicates the position of the 
Fermi level (Ef = eV). (a) Only Vo are present in the cell, 
(b) V03 vacancy co-doped with N. Note that in this latter 
case the system is magnetic. 



that the system is non-magnetic and remains effectively 
semi-conducting. If we assume for the moment that the 
clusters of Hf atoms co-ordinating the V04 and the V03 
sites form a perfect tetrahedron and a perfect trigonal 
plane respectively, resulting in local Td and C^h symme- 
tries, then the single particle molecular orbital ground 
state in each case would be a completely symmetric and 
non-degenerate al and aV singlet respectively. Similarly 
the higher lying cr-bonding single particle states would 
be a t% triplet in the V04 case and a e' doublet in the 
V03 case. Deviation from perfect symmetry, as in the 
actual case, of course means that the degeneracies of the 
excited states are lifted and energy levels suitably re- 
ordered. Nevertheless, the ground state remains an or- 
bital singlet. The higher lying states are well separated 
in energy from the low-lying singlet level. Configurations 
with one electron promoted to the higher lying states lie 
prohibitively higher in energy relative to the ground state 
thus ruling out electron promotion. Since two electrons 
occupying an orbital singlet anti-align their spins, the 
resulting ground state is non-magnetic. 

We then investigate the effect of partial p-doping of the 
system, by substituting an electron acceptor like N at an 
O site. In figure [2h we present the DOS for V03 only, 
bearing in mind that the situation for V04 is essentially 
the same. In this case one electron is removed from the 
Vb level, which therefore spin splits with a magnetic mo- 
ment of 1 /is per vacancy. We then check whether these 
localized moments interact with each other, by doubling 
our supercell and comparing the total energies of the fer- 
romagnetic and anti-ferromagnetic alignment of the mag- 
netic moments on two seperate vacancy sites. We find 
no difference in the total energy and so we conclude that 
the isolated moments are not coupled. The system can at 
best be paramagnetic and thus the Vo themselves cannot 
possibly be behind the observed ferromagnetism. 

We now move on to investigate the hypothesis of ferro- 
magnetic order due to Vm- In monoclinic Hf02 (badde- 
leyite), each Hf atom is co-ordinated by seven O atoms 
(see figure Of the seven O atoms, three are of one 
type which we label '01' and the remaining four of a 
second type which we label '02'. The two types of O 
atoms differ in their Hf co-ordination number in the crys- 
tal. We further label the remaining O atoms not directly 
co-ordinating the vacancy site as type '03'. In the per- 
fect crystal, the O 2p levels are fully filled and form the 
bulk of the valence band. Since Hf is a cation with a 4+ 
valence, the removal a neutral Hf atom introduces four 
empty states among the oxygen 2p levels. The spin oc- 
cupation of these four states established whether or not 
the system is magnetic. 

At this point we present our LSDA results for the 
ground state of Hf deficient Hf02- We consider a 2x2x2 
supercell containing 96 atoms. Only one Hf vacancy is 
introduced in the supercell and upon relaxation the O 
atoms are seen to move outwards around the vacancy site 
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by about 0.15 A. In figure. |3|we present the DOS for the 
fully relaxed case. The valence band is clearly spin split, 
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FIG. 3: Spin resolved DOS for one Vm is a 96 atoms HfCh 
supercell. Ef is set at 0.0 eV. 




FIG. 5: Localization plot for the magnetic moment inside the 
supercell before and after relaxation. The quantity presented 
is pol(R) = fft +dR dr JdQ[pi(r) — Pi(r)] against R. pol(R) 
is the spin polarization from a shell of thickness dR at radius 
R where R is measured radially outward from the vacancy 
site. Note the redistribution of the magnetic moment towards 
outer shells upon relaxation. 



with the compensating holes mostly confined to the down 
spin states resulting in an almost half metallic ground 
state with a magnetic moment of 3.52 /zb per vacancy. 
Prior to relaxation the system is completely half metallic 
with an integer magnetic moment of 4 fi^ per vacancy. 
The relaxation involves a considerable redistribution of 
the hole density over the O atoms around the vacancy site 
as shown if figure 0] with the magnetism coming predom- 
inantly from the 01 type atoms. The observed charge re- 
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FIG. 4: Isosurface of the local density of states for the high- 
est, minority-spin defect level of Vm- After relaxation the 
charge density (hole density) clumps mainly on the three Ol 
type O atoms. As a result the 02 type atoms lose most of 
their spin polarization. The central atom is a dummy atom 
included only for clarity. 

ordering is driven by large scale re-hybridization, upon 
relaxation, of the orbitals constituting the impurity lev- 
els with the crystalline surroundings. We find that the 
holes clump together on and around the Ol atoms and 
their nearest 03 type neighbours. Figure shows the 
localization of the magnetic moment around the vacancy 
site, where this outward spread in the polarization due 
to relaxation is evident. 

Having studied the real space distribution of the mag- 
netic moment, we now look at the symmetry aspect of 
the defect states. Considering the seven O atoms co- 
ordinating the vacancy site as a single molecular clus- 



ter, we find a trivial C\ point group. Thus the sin- 
gle particle molecular orbitals generated from combin- 
ing the O atomic orbitals form a set of non-degenerate 
levels, and the high spin state arises from the singly oc- 
cupied four topmost molecular orbitals (see insets of fig- 
ure EJ . Interestingly the calculated magnetic moment 
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FIG. 6: Magnetic moment per vacancy as a function of 
the average distance from the vacancy site of the seven co- 
ordinating Oxygen atoms. In the two insets schematic repre- 
sentations of the high spin S=2, and non-magnetic S=0 state. 

(^4^b) cannot originate from a single orbitally degen- 
erate molecular level even in the most symmetric octa- 
hedral molecule. In fact the largest orbital degeneracy 
allowed is just three fold, and a total spin S=2 for a con- 
figuration of four spin 1/2 particles is ruled out by the 
Pauli principle. This means that in any case a set of non 
degenerate orbitals must be involved in the high spin con- 
figuration. As non-degenerate molecular orbitals differ in 
single particle properties (kinetic and ionic energies), the 
final configuration of the four spin 1/2 particles is de- 
cided not just by mutual Coulomb repulsion or exchange 
energy but by an interplay between all contributions in 
the Hamiltonian. 

Looking at the relative energetics of different possible 
electronic configurations of the defect levels helps to see 
what stabilizes the high spin ground state. To this end, 
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wc perform fixed spin moment calculations, for the two 
limiting cases of S=2 and S=0. From the calculations 
it emerges that the high spin configuration is expensive 
with regards to kinetic and Hartree energies but the gain 
in exchange and ionic energy is enough to stabilize it. 
It is likely that the higher lying molecular orbitals are 
strongly anti-bonding in character with the electrons in 
these levels being localized more on the ions. This results 
in higher kinetic energies and lower ionic potential ener- 
gies. In fact, one expects that analogously to what hap- 
pens in molecules, if the O atoms around the vacancy are 
artificially squeezed in towards the vacancy, thus driving 
the system highly kinetic, the higher lying defect levels 
would be emptied out accompanied by a fall in the mag- 
netic moment. This is indeed seen to be the case as shown 
in figure EJ As it turns out, at the equilibrium distance 
in the crystal (~2.28A) , the magnetic moment is close 

to ifiB- 

Next we address the question of the magnetic coupling 
between vacancies, by calculating the total energy of a 
supercell containing two Vkf and comparing the energy 
for the ferromagnetic (Efm) an d antiferromagnetic state 
(Eaf) (see Table HJ. Clearly the ferromagnetic align- 



Vm concentration 


AE (meV) 


2.08% (1) 


-113.66 


1.38% (2) 


-50.56 


1.04% (3) 


-13.81 



TABLE I: Energy differences AE = Efm — Eaf between fer- 
romagnetic and antiferromagnetic alignment of the magnetic 
moments of two different vacancy sites. The first column in- 
dicates the defect concentration and the number in bracket 
the number of Hf sites separating the two vacancies. 



ment is always energetically favorable, and most remark- 
ably the coupling appears rather strong, in particular for 
short Vm — Vm distances. This leads us to attribute the 
observed ferromagnetism in HfC-2 thin-films |l| to Hf va- 
cancies. Moreover the large values of AE suggest Curie 
temperatures above room temperature at large enough 
concentrations . 

In order to have a better understanding of the ori- 
gin of this ferromagnetic ordering we have performed an 
extensive study of the charge distribution by means of 
Miilliken population analysis. The main features are: (a) 
the O atoms in the cell are polarized to different degrees 
depending on their orientation and distance relative to 
the vacancy site but always with the same sign, (b) the 
Hf atoms in the cell are also polarized but importantly, 
the sign of polarization is opposite to that of the O, and 
(c) the total polarization of all the O atoms in the cell 
is 3.92/zb and that for the Hf atoms is -0.40/iB, leaving 
a moment of the cell of 3.52/iB- This suggests that the 
magnetic coupling between the O atoms in the cell is me- 
diated by minority spin electron derealization across the 



Hf bridge connecting the O atoms. This applies also to 
O atoms belonging to two different Vm sites. The der- 
ealization is larger when local moments on the two Vm 
are ferromagnetically aligned resulting in lower Kinetic 
and Exchange energies relative to the anti-ferromagnetic 
case. 

In summary we have performed DFT calculations 
investigating the possibility of intrinsic defect driven 
ferromagnetism in Hf02- Oxygen vacancies form 
non-magnetic impurity levels unless p-type co-dopants 
are present. However, in this case the magnetic in- 
teraction is negligible, ruling out the hypothesis of d° 
ferromagnetism Q. In contrast Hf vacancies show a 
high spin state with an associated magnetic moment 
of ^3.5 fiB- These are ferromagnetically coupled via 
minority spin electron derealization across the bridging 
Hf sites, with a large coupling strength suggesting high 
Curie temperatures. Thus a new direction for intrinsic 
defect ferromagnetism is evident based on the following 
facts: 1) symmetry driven orbital degeneracy is not 
a pre-requisite for the existence of a high-spin defect 
ground state, 2) a set of closely spaced single particle 
levels together with strong exchange might be sufficient 
for the same, 3) suitable mixing of defect states with the 
crystalline environment can lead to ferromagnetic inter- 
defect coupling. These findings suggest that a wider 
class of systems, not restricted by symmetry and free 
from the possibility of Jahn- Teller like distortions, might 
actually be open to intrinsic defect ferromagnetism. 

We thank Prof. J.M.D Coey for useful discussions. 
This work is funded by the Science Foundation of Ireland 
under the grant SFI02/LN1/I175. 



[1] M. Venkatesan, C. B. Fitzgerald, J. M. D. Coey, Nature 

430, 630 (2004) Brief Communications. 
[2] N.I. Medvedeva et.al, Phys. stat. sol. (b) 160, 517 (1990). 
[3] X. Zhao and D. Vanderbilt, Phys. Rev. B 65 233106 

(2002). 

[4] A. S. Foster, F. Lopez Gejo, A. L. Shluger and R. M. 

Nieminen, Phys. Rev. B 65, 174117. 
[5] W.L. Scopel, A.J.R. da Silva, W. Orellana and A. Fazzio, 

Appl. Phys. Lett. 84, 1492 (2004). 
[6] R. Monnier and B. Delley, Phys. Rev. Lett. 87 157204 

(2001). 

[7] I. S. Elfimov, S. Yunoki, and G. A. Sawatzky, Phys. Rev. 

Lett. 89, 216403 (2002). 
[8] A. Zywietz, J. Furthmiiller, and F. Bechstedt, Phys. Rev. 

B 62, 6854 (2000). 
[9] D.M. Ceperley and B.J. Alder, Phys. Rev. Lett. 45, 566 

(1980) 

[10] J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Jun- 
quera, P. Ordejon and D. Sanchez-Portal, J. Phys. Cond. 
Matter 14, 2745 (2002). 



